Tailoring the facet distribution on copper with chloride.

Electrocatalytic reactions are sensitive to the catalyst surface structure. Therefore, finding methods to determine active surface sites with different geometry is essential to address the structure-electrocatalytic performance relationships. In this work, we propose a simple methodology to tune and quantify the surface structure on copper catalysts. We tailor the distribution and ratio of facets on copper by electrochemically oxidizing and reducing the surface in chloride-rich aqueous solutions. We then address the formation of new facets with voltammetric lead (Pb) underpotential deposition (UPD). We first record the voltammetric lead UPD on different single facets, which have intense peaks at different potential values. We use this data to decouple each facet peak-contribution in the lead (Pb) UPD curves of the tailored and multifaceted copper surfaces and determine the geometry of the active sites. We combine experiments with density functional theory (DFT) calculations to assess the ligand effect of chloride anions on the copper facet distribution during the surface oxidation/electrodeposition treatment. Our experiments and Wulff constructions suggest that chloride preferentially adsorbs on the (310) facet, reducing the number of (111) sites and inducing the growth of (310) or n(100) × (110) domains. Our work provides a tool to correlate active sites with copper geometries, which is needed to assess the structure-performance relationships in electrocatalysis. We also demonstrate an easy method for selectively tailoring the facet distribution of copper, which is essential to design a well-defined nanostructured catalyst.

Tuning the chemical composition of binary alloy nanoparticles to prevent their dissolution.

The dissolution of nanoparticles under corrosive environments represents one of the main issues in electrochemical processes. Here, a model for alloying and protecting nanoparticles from corrosion with an anti-corrosive element (e.g. Au) is proposed based on the hypothesis that under-coordinated atoms are the first atoms to dissolve. The model considers the dissolution of atoms with coordination number ≤6 on A-B nanoparticles with different sizes, shapes, chemical compositions, and exposed crystallographic orientations. The results revealed that the nanoparticle's size and chemical composition play a key role in the dissolution, suggesting that a certain composition of an element with corrosive resistance could be used to protect nanoparticles. DFT simulations were performed to support our model on the dissolution of four types of atoms commonly found on the surface of Au0.20Pd0.80 binary alloys - terrace, edge, kink, and ad atoms. The simulations suggest that the less coordinated ad and kink Pd atoms on Au0.20Pd0.80 alloys are dissolved in a potential window between 0.26-0.56 V, while the rest of the Pd and Au atoms are protected. Furthermore, to show that a corrosion-resistant element can indeed protect nanoparticles, we experimentally investigated the electrochemical dissolution of immobilized Pd, Au0.20Pd0.80, and Au0.40Pd0.60 nanoparticles in a harsh environment. In line with the dissolution model, the experimental results show that an Au molar fraction of the nanoparticle of 0.20, i.e., Au0.20Pd0.80 binary alloy, is a good compromise between maximizing the active surface area (Pd atoms) and corrosion protection by the inactive Au.

Structure and energetics of liquid water-hydroxyl layers on Pt(111).

The interactions between liquid water and hydroxyl species on Pt(111) surfaces have been intensely investigated due to their importance to fuel cell electrocatalysis. Here we present a molecular dynamics study of their structure and energetics using an ensemble of neural network potentials, which allow us to obtain unprecedented statistical sampling. We first study the energetics of hydroxyl formation, where we find a near-linear adsorption energy profile, which exhibits a soft and gradual increase in the differential adsorption energy at high hydroxyl coverages. This is strikingly different from the predictions of the conventional bilayer model, which displays a kink at 1/3ML OH coverage indicating a sizeable jump in differential adsorption energy, but within the statistical uncertainty of previously reported ab initio molecular dynamics studies. We then analyze the structure of the interface, where we provide evidence for the water-OH/Pt(111) interface being hydrophobic at high hydroxyl coverages. We furthermore explain the observed adsorption energetics by analyzing the hydrogen bonding in the water-hydroxyl adlayers, where we argue that the increase in differential adsorption energy at high OH coverage can be explained by a reduction in the number of hydrogen bonds from the adsorbed water molecules to the hydroxyls.

Addressing Dynamics at Catalytic Heterogeneous Interfaces with DFT-MD: Anomalous Temperature Distributions from Commonly Used Thermostats.

Density functional theory-based molecular dynamics (DFT-MD) has been widely used for studying the chemistry of heterogeneous interfacial systems under operational conditions. We report frequently overlooked errors in thermostated or constant-temperature DFT-MD simulations applied to study (electro)catalytic chemistry. Our results demonstrate that commonly used thermostats such as Nosé-Hoover, Berendsen, and simple velocity-rescaling methods fail to provide a reliable temperature description for systems considered. Instead, nonconstant temperatures and large temperature gradients within the different parts of the system are observed. The errors are not a "feature" of any particular code but are present in several ab initio molecular dynamics implementations. This uneven temperature distribution, due to inadequate thermostatting, is well-known in the classical MD community, where it is ascribed to the failure in kinetic energy equipartition among different degrees of freedom in heterogeneous systems (Harvey et al. J. Comput. Chem. 1998, 726-740) and termed the flying ice cube effect. We provide tantamount evidence that interfacial systems are susceptible to substantial flying ice cube effects and demonstrate that the traditional Nosé-Hoover and Berendsen thermostats should be applied with care when simulating, for example, catalytic properties or structures of solvated interfaces and supported clusters. We conclude that the flying ice cube effect in these systems can be conveniently avoided using Langevin dynamics.

How to extract adsorption energies, adsorbate-adsorbate interaction parameters and saturation coverages from temperature programmed desorption experiments.

We present a scheme to extract the adsorption energy, adsorbate interaction parameter and the saturation coverage from temperature programmed desorption (TPD) experiments. We propose that the coverage dependent adsorption energy can be fit using a functional form including the configurational entropy and linear adsorbate-adsorbate interaction terms. As one example of this scheme, we analyze TPD of CO desorption on Au(211) and Au(310) surfaces. We determine that under atmospheric CO pressure, the steps of both facets adsorb between 0.4-0.9 ML coverage of CO*. We compare this result against energies obtained from five density functionals, RPBE, PBE, PBE-D3, RPBE-D3 and BEEF-vdW. We find that the energies and equilibrium coverages from RPBE-D3 and PBE are closest to the values determined from the TPD.

Solvation at metal/water interfaces: An ab initio molecular dynamics benchmark of common computational approaches.

Determining the influence of the solvent on electrochemical reaction energetics is a central challenge in our understanding of electrochemical interfaces. To date, it is unclear how well existing methods predict solvation energies at solid/liquid interfaces, since they cannot be assessed experimentally. Ab initio molecular dynamics (AIMD) simulations present a physically highly accurate, but also a very costly approach. In this work, we employ extensive AIMD simulations to benchmark solvation at charge-neutral metal/water interfaces against commonly applied continuum solvent models. We consider a variety of adsorbates including *CO, *CHO, *COH, *OCCHO, *OH, and *OOH on Cu, Au, and Pt facets solvated by water. The surfaces and adsorbates considered are relevant, among other reactions, to electrochemical CO2 reduction and the oxygen redox reactions. We determine directional hydrogen bonds and steric water competition to be critical for a correct description of solvation at the metal/water interfaces. As a consequence, we find that the most frequently applied continuum solvation methods, which do not yet capture these properties, do not presently provide more accurate energetics over simulations in vacuum. We find most of the computed benchmark solvation energies to linearly scale with hydrogen bonding or competitive water adsorption, which strongly differ across surfaces. Thus, we determine solvation energies of adsorbates to be non-transferable between metal surfaces, in contrast to standard practice.

Fingerprint Voltammograms of Copper Single Crystals under Alkaline Conditions: A Fundamental Mechanistic Analysis.

A critical step toward the systematic development of electrocatalysts is the determination of the microscopic structure and processes at the electrified solid/electrolyte interface. The major challenges toward this end for experiment and computations are achieving sufficient cleanliness and modeling the complexity of electrochemical systems, respectively. In this sense, benchmarks of well-defined model systems are sparse. This work presents a rigorous joint experimental-theoretical study on the single-crystal (SC) Cu/aqueous interface. Within typical computational uncertainties, we find quantitative agreement between simulated and experimentally measured voltammograms, which allows us to unequivocally identify the *OH adsorption feature in the fingerprint region of Cu(110), Cu(100), and Cu(111) SCs under alkaline conditions. We find the inclusion of hydrogen evolution reaction kinetics in the theoretical model to be crucial for an accurate steady-state description that gives rise to a negligible H* coverage. A purely thermodynamic description of the H* coverage through a Pourbaix analysis would incorrectly lead to a H* adsorption peak. The presented results establish a fundamental benchmark for all electrochemical applications of Cu.

Energy-entropy competition in cation-hydroxyl interactions at the liquid water-Pt(111) interface.

Electrochemical reaction rates are sensitive to interactions between electrolyte cations and adsorbed reaction intermediates, e.g., cation-*OH interactions in the oxygen reduction reaction on platinum. Here, we calculate the free energy interaction between adsorbed *OH and K+/Li+ situated at the liquid water-Pt(111) interface using ab initio molecular dynamics (AIMD) and metadynamics. Li+ stabilizes *OH by 0.1 ± 0.1 eV and K+ destabilizes *OH by 0.1 ± 0.1 eV, in qualitative agreement with experimental cyclic voltammogram (CV) measurements. In contrast, the internal energy of *OH is stabilized by 0.3 eV and 0.4 eV for Li+ and K+, respectively. This demonstrates that entropy significantly destabilizes cation-*OH interactions and is vital in order to understand even the relative influence of cations at interfaces.

OH formation and H2 adsorption at the liquid water-Pt(111) interface.

The liquid water-Pt(111) interface is studied with constant temperature ab initio molecular dynamics to explore the importance of liquid water dynamics of catalytic reactions such as the oxygen reduction reaction in PEM fuel cells. The structure and energetics of hydroxyls formed at the liquid water-Pt(111) interface are found to be significantly different from those of the hydroxyl formed on a bare Pt(111) surface and the hydroxyl formed on a Pt(111) surface with a static water layer. We identify 1/12 ML *OH, 5/12 ML *OH and 2/3 ML *OH as particularly stable hydroxyl coverages in highly dynamic liquid water environments, which - contrary to static water-hydroxyl models - contain adjacent uncovered Pt sites. Atomic surface oxygen is found to be unstable in the presence of liquid water, in contrast to static atomic level simulations. These results give an improved understanding of hydroxide and surface oxide formation from Pt(111) cyclic voltammetry and allow us to draw detailed connections between the electrostatic potential and the interface structure. The study of hydrogen adsorption at the liquid water-Pt(111) interface finds competitive adsorption between the adsorbed hydrogen atoms and water molecules. This does not adhere with experimental observations, and this indicates that the Pt(111) surface has to be negatively charged for a correct description of the liquid water-Pt(111) interface at potentials where hydrogen adsorption occurs.

Catechol and HCl Adsorption on TiO2(110) in Vacuum and at the Water-TiO2 Interface.

Coadsorbed water is often unavoidable in electrochemistry and low-temperature catalysis. In addition, water influences the adsorption of biomolecules on surfaces. We use ab initio DFT molecular dynamics and ground-state calculations to study the adsorption of HCl and catechol on the rutile TiO2(110) surface and at a water-rutile interface. We find that a coadsorbed water film reduces the adsorption energy of both catechol and HCl significantly because water molecules must be displaced from the surface before catechol or HCl can adsorb. The adsorption energy of catechol (or HCl) at the water-rutile interface can be estimated as the adsorption energy in vacuum minus the energy to remove two water molecules (respectively, one water molecule) from the rutile surface in vacuum and place them in liquid water. This estimate predicts the effect of a surface water film on adsorption without the need of molecular dynamics.

Density functional theory and chromium: Insights from the dimers.

The binding in small Cr clusters is re-investigated, where the correct description of the dimer in three charge states is used as criterion to assign the most suitable density functional theory approximation. The difficulty in chromium arises from the subtle interplay between energy gain from hybridization and energetic cost due to exchange between s and d based molecular orbitals. Variations in published bond lengths and binding energies are shown to arise from insufficient numerical representation of electron density and Kohn-Sham wave-functions. The best functional performance is found for gradient corrected (GGA) functionals and meta-GGAs, where we find severe differences between functionals from the same family due to the importance of exchange. Only the "best fit" from Bayesian error estimation is able to predict the correct energetics for all three charge states unambiguously. With this knowledge, we predict small bond-lengths to be exclusively present in Cr2 and Cr2 (-). Already for the dimer cation, solely long bond-lengths appear, similar to what is found in the trimer and in chromium bulk.

DFT+U study of defects in bulk rutile TiO(2).

We present a systematic study of electronic gap states in defected titania using our implementation of the Hubbard-U approximation in the grid-based projector-augmented wave density functional theory code, GPAW. The defects considered are Ti interstitials, O vacancies, and H dopants in the rutile phase of bulk titanium dioxide. We find that by applying a sufficiently large value for the Hubbard-U parameter of the Ti 3d states, the excess electrons localize spatially at the Ti sites and appear as states in the band gap. At U=2.5 eV, the position in energy of these gap states are in fair agreement with the experimental observations. In calculations with several excess electrons and U=2.5 eV, all of these end up in gap states that are spatially localized around specific Ti atoms, thus effectively creating one Ti(3+) ion per excess electron. An important result of this investigation is that regardless of which structural defect is the origin of the gap states, at U=2.5 eV, these states are found to have their mean energies within a few hundredths of an eV from 0.94 eV below the conduction band minimum.